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1. INTRODUCTION 

Currently, there are different methods and algorithms for image processing that are increasingly 
used in the medical field to solve various problems: search for anomalies in medical images, classification of 
X-rays, and assistance in the analysis. Finding the key points is an important issue for solving the tasks in 
various fields, such as video surveillance, parameterization, identification of objects, and image stitching. 
Key points should be: i) uniquely identifiable in some neighborhoods and image as a whole; ii) invariant to 
affine transformations; iii) stable relative to noise; iv) effective for parameterization and identification of 
objects in the set of analyzed images. 

As a key points for images can be: pixel, line, and intersection of lines—corners. There are a number 
of methods for finding corners in an image: Harris, [1], Shi-Tomasi [2], and FAST [3]. Corner detectors 
most often work with the brightness component of the pixel value. To the input of such a detector, working 
with a luminance component, black and white image served, as a result, and the output is formed by a list of 
possible angles with the degree of similarity. 

Points with a degree of similarity greater than the threshold are defined as angles and less than the 
threshold—discarded. The disadvantage of using such detectors in the above tasks is that: that they detect the 
area, where the corner point is found, but do not always locate its coordinates, such as the corner point as a 
result of the work of these algorithms may be outside the line. Therefore, after application of these detectors 
is necessary to specify the location of each angle, for example using a mask. To eliminate these drawbacks, it 
is necessary to use angles and lines as key points. 
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2. METHOD 

For tasks such as dividing paths into separate straight lines and stitching images, corners are used as 
key points. Corner can be classified as shown in: i) L-form (with different angles); 11) Y-form (with different 
angles) or usually like T-form; iii) X-form (with different angles). As shown in Figure 1 developed an 
algorithm for corner detection in aerospace image by using an expandable mask. The proposed algorithm is 
based on the sequential application of a mask at each pixel of a selected contour, not a straight line, and 
except for endpoints of the contour. In case of an ambiguous solution for a pixel, it is additionally checked by 
expanding the mask to the size pixels. 


L Y TTX 


Figure 1. Types of corners in the image 


The input of the algorithm is an image 7 = ||i(y, lo = OY=1, x = OX=1) where i(y, x) = 0.255- 
the brightness value of a pixel in the image, Y, X — the dimensions of the image vertically and horizontally. 
The algorithm is consisting of the following steps. 

— Formation of a matrix of binary images Iz = |lig(y, Dll = 07I, x -0X T} for the input image / using 
the canny contour filtering algorithm [4], where ipg(y, x) =1 for a pixel belonging to the path, 
ig(y, x) = 0 for the background pixel. 

— Segmentation of contour lines. Segmentation is carried out by the method (region growing (RG)) 
[5]-[8]. Each contour pixel ig(y, x) = 1 the contour number is assigned, to which it belongs. As a 
result, a matrix of contours is formed S = ||s(7)||,,-z77, and a matrix of the number of endpoints in each 
contour K = ||k(m)||,-¢n- Where s(n)- coordinates of contour pixels, represented as matrices 
X(n) = ||x(n, M!l(c=0e=1)> Y(n) = |ly(n, Nl(c=0e=1)> k(n) is number of the endpoint, N is number of 
found contours, C is number of pixels in the contour shows in Figures 2 (a) and (b). 


(a) (b) 


Figure 2. Segmentation of contour lines (a) test image and (b) segmentation (RG) for test image 


— Normalization of selected contours by thickness. Contour normalization is carried out using the method 
of contour line thickness normalization [9]-[11]. In the process of normalization X(n), Y (n) pixels are 
removed from the contour coordinate matrices, which visually and physically make the line thicker. As 
a result, s(n), outlines are formed with a thickness of one pixel. 

— Contour analysis. The endpoints of the contours are analyzed first k(n). If k(n) = 2- the decision is 
made that the contour s(n) is a line and the form factor is calculated f to check if the contour is curved 
s(n) [12]-[15]. If f = [0.8,1.2]— the decision is made that the contour s(n) is straight lines and 
therefore has no corners. Therefore, this contour is removed from the matrix S. 

— Search corner points. As a result of the execution of the algorithm, angular matrices are formed X, = 
Ix(¢, MIlj-a2- Y, = ly Mly- where x(g,n), y(g,n)— the coordinates of the corner points, g 
issequence number of the coordinates, G is the number of corner points found, n is contour number to 
which the corner point belongs [16]-[19]. 
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a. Formatting mask M; = lm, in Figure 3, where j = 1,16 is sequence number of the 


Dll y- 


matrix, 
1, for j = {1,2,9}, 1 for j = {5,6,13,14}, 
(0,0) = Se Oy a hp 
EAER) k for j = (3.8, 10,16}, 0) 10, for j = {1.4,7. 12, 15,16}, 
1, for j = {1,4,11,16}, 1, for j = {6,7,11,12}, 
0, for j = {2,3,5..10, 125-15, }, 0, for j = {1..5,8..10,13..16}, 
m;(1,1)=1, 
1, for j = {5,8..10}, 1, for j = {2,3,10,14}, 
m;(1,2) -| — „= "y0) = P TO Ties 
i 0, for j= {1..4,6,7,11..16}, ” 0, for j = {1,4..9,11..13, 15,16}, 
1, for j = {7,8,15,16}, 1, for j = {3,4,12,13}, 
0, for j = {1..6,9..14}, 0, for j = {1,2,5..11,14..16}. 
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Figure 3. Masks for searching corner points in contours 


b. Forming corner matrices Y}, X4. The neighborhood of each pixel of the contour s(n) with coordinates 
y(n,c),x(n,c) matrix Y(n),X(n) is checked for compliance with one of 16 corner masks Mj. If the 
neighborhood of a pixel with coordinates y(n,c),x(n,c) matches masks M,..Mg — go to the step d. If 
the neighborhood of the pixel with coordinates y(n, c), x(n, c) corresponds to the masks Mo. . M16- go 
to step c [20]-[23]. 

c. Clarification. For a contour pixel with coordinates y(n,c),x(n,c), the neighborhood of which 
corresponds to matrices My..M,,, refinement is performed by expanding the matrix to the size 5 x 5 
pixels as shown in Figure 4. If the neighborhood of the pixel with coordinates y(n,c),x(n,c) 
corresponds to the expanded matrix (the presence of at least one contour point on each side is checked 
in the area indicated in gray in Figure 2) go to step d. 


13 14 15 16 


Figure 4. Masks for searching corner points in contours 
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d. Pixel coordinates y(n, c), x(n, c), defined as a corner point are entered into matrices Y,,X, accordingly 
[24]-[26]. As a result of the algorithm, coordinate matrices are formed Y,, X, containing the coordinates 
of the corner points Figure 5. 
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Figure 5. Corner detection in areospace image by using an expandable mask 


3. RESULTS AND DISCUSSION 

The developed algorithm is implemented in the C** programming language using the open computer 
vision (CV) library. For comparative evaluation algorithm implemented most known methods for angles 
detection—Harris. The experiment was tested on a computer with the following technical characteristics: 
processor—intel (R) core (TM) 15-2320 central processing unit (CPU) 3.0 GHz; random acces memory-4 
(RAM-4) GB; system type—64-bit operating system, x64 processor; and operating system—windows 7. 

To evaluate the operation of the algorithms, grayscale images were divided into 5 classes according 
to the type of the brightness of the image histogram. One image from each class was taken for testing. The 
test image and their histograms are shown in Figures 6(a) to (e). 
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(e) 
Figure 6. Test images and their histograms of different classes (a) —1, (b) —2, (c) —3, (d) —4, and (e) —5 


For quantitative evaluation of stability S(x) the number of detected corners was compared, retained 


their location after changing in the brightness, contrast, and in the image, rotation as shown in Figures 7(a) to 
(1), Figures 8(a) to (1), Figures 9(a) to (m), and u = ||b, c, a|| with the number of detected corners in the test 
image. To evaluate the stability of the corner detection, their number C;() on the image i, when one of its 
parameters was changed, it was normalized relative to the test image according to: 


sa =e (1) 


where C, is the number of detected corners in the test image. u = ||b,c,a|], b = [—50,50],%, 
c = [-50,50], %, a = [—90, 90], ° 
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(i) Q) (k) (D 


Figure 7. Changing in the brightness (a) test image, (b) 0, (c) 10, (d) 20, (e) 30, (£) 40, (g) 50, (h) -10, (i) -20, 
(j) -30, (k) -40, and (1) -50 


Figure 8. Changing in the contrast (a) test image, (b) 0, (c) 10, (d) 20, (e) 30, (£) 40, (g) 50, (h) -10, (i) -20, 
(j) -30, (k) -40, and (1) -50 
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Figure 9. Changing in the rotation, (a) test image, (b) 15, (c) 30, (d) 45, (e) 60, (f) 75, (g) 90, (h) -15, (i) -30, 
(j) -45, (k) -60, (1) -75, and (m) -90 


As shown in Figures 10(a) to (0), the developed algorithm wins up to 44% over the Harris method 
when changing the brightness from -50 to 50 for the 5 class of the test images Figure 10(a) and the proposed 
algorithm is slower than Harris method when the brightness changes from -50 to -10 for the 3rd class of 
images up to 60% Figure 10(b). With a decrease in contrast in the image, the proposed algorithm is slower 
than Harris method by up to 90%, and with an increase—it wins up to 60% Figure 10(c). When rotating the 
image, the developed algorithm is more stable for the 5" class of images Figure 10(d) and less stable for 2 
and 4 classes of images Figures 10(e), (f). 
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Figure 10. Graph for evaluating the stability of corner detection algorithm for different classes and changing 
parameters of the image (a) class 1, change in the brightness, (b) class 1, change in the contrast, (c) class 1, 
rotation, (d) class 2, change in the brightness, (e) class 2, change in the contrast, (f) class 2, rotation, (g) class 
3, change in the brightness, (h) class 3, change in the contrast, (i) class 3, rotation, (j) class 4, change in the 
brightness, (k) class 4, change in the contrast, (1) class 4, rotation, (m) class 5, change in the brightness, (n) 

class 5, change in the contrast, and (o) class 5, rotation 
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The average running time for each of the test images is shown in Table 1. From the Table 1. It is 


seen that the average running time of the presented method is 2.5—13.3 times less than the average running 
time of the Harris method. 


4. 


Table 1. Average running time of the test algorithms (s) 


Algorithms Image 1 _Image2 _Image3 Image 4_ Image 5 
Proposed algorithm 1.3 1.2 1.3 0.5 0.6 
Harris method 3.3 16 9.8 2.4 2.4 


CONCLUSION 
An algorithm for corner detection in aerospace image based on an expanding mask has been 


developed. It is shown that the proposed algorithm is up to 60% more stable than the Harris method with 
increasing contrast. The proposed algorithm is more stable than the Harris method when rotating images of 
the 5 class, but inferior to Harris method for all other classes of images when rotating at different angles. The 
disadvantages of the proposed algorithm also include the loss of stability with a decrease in contrast. The 
above disadvantages of the proposed algorithm are compensated for by the high performance relative to the 
Harris method. It is shown that the proposed algorithm faster Harris's method by 13.3 times in time. 
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